Rfast (version 1.5)

Fitting a beta distribution via Newton-Rapshon: Fitting a beta distribution via Newton-Rapshon

Description

Fitting a Dirichlet distribution via Newton-Rapshon.

Usage

beta.mle(x)

Arguments

x

A numerical vector with proportions, i.e. numbers in (0, 1) (zeros and ones are not allowed).

Value

A list including: A list including:

Details

Maximum likelihood estimation of the parameters of a beta distribution is performed via Newton-Raphson. The x values and 1 - x become a two-column matrix and the function diri.nr2 is implemented.

See Also

diri.nr2

Examples

x <- rbeta(1000, 1,4)
system.time( for(i in 1:1000) beta.mle(x) )